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Calibrating [Oil] star- format ion rates at z>l from dual 
Ha- [Oil] imaging from HiZELS 



(N 

o 

(N 
O 

Q 

o 

(N 



o 

u 

a,: 
6 

CO 

o3 



> 

in 
o 

On 
(N 



IS 



Masao Hayashi 1 *, David Sobral 2 , Philip N. Best 3 , Ian Smail 4 , 
and Tadayuki Kodama 1 ' 5 

1 Optical and Infrared Astronomy Division, National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan 

2 Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, THe Netherlands 

^SUPA, Institute for Astronomy, Royal Observatory of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK 

4 Institute for Computational Cosmology, Durham University, South Road, Durham, DH1 3LE, UK 

5 Subaru Telescope, National Astronomical Observatory of Japan, 650 North A'ohoku Place, Hilo, HI 96720, USA 



21 December 2012 



ABSTRACT 

We investigate the relationship between Ha and [On](A3727) emission in faint star- 
forming galaxies at z = 1.47 with dust uncorrected star formation rates (SFRs) down 
to 1.4 M Q yr _1 , using data in two narrow-bands from WFCAM/UKIRT and Suprime- 
Cam/Subaru. A stacking analysis allows us to investigate Ha emission flux from 
bright [O n] emitters as well as faint ones for which H a is not individually detected, 
and to compare them with a large sample of local galaxies. We find that there is a 
clear, positive correlation between the average H a and [O n] luminosities for [O n] 
emitters at z = 1.47, with its slope being consistent with the local relation. [On] 
emitters at z = 1.47 have lower mean observed ratios of Ha/ [On] suggesting a small 
but systematic offset (at 2.8a significance) towards lower values of dust attenuation, 
^4r"c(~ 0.35, than local galaxies. This confirms that [On] selection tends to pick up 
galaxies which are significantly less dusty on average than Ha selected ones, with the 
difference being higher at z = 1.47 than at z = 0. The discrepancy of the observed line 
ratios between [On] emitters at z = 1.47 and the local galaxies may in part be due to 
the samples having different metallicities. However, we demonstrate that metallicity 
is unlikely to be the main cause. Therefore, it is important to take into account that 
the relations for the dust correction which are derived using H a emitter samples, 
and frequently used in many studies of high-z galaxies, may overestimate the intrinsic 
SFRs of [O n] -selected galaxies, and that surveys of [On] emission galaxies are likely 
to miss dusty populations. 

Key words: galaxies: high-redshift - galaxies: evolution. 



1 INTRODUCTION 

Star formation rate (SFR) is one of the most important 
properties to characterise the growth of a galaxy, and the 
star formation history of the Universe provides us with a 
fundamental insight into galaxy evolution. Several surveys 
have been conducted to reveal the star formation activity 
in the distant Universe by making use of various indicators 
such as ultraviolet (UV) luminosity, nebular emissions such 
as Ha and [O Il](A3727), and infrared (IR) radiation. These 
have been used to estimate SFR of galaxies over a wide range 
of redshifts; e.g., UV: Ouchi et al. (2009); Bouwens et al. 
(2011), nebular lines: Shioya et al. (2008); Ly et al. (2007, 



E-mail: masao.hayashi@nao.ac.jp 



2011); Geach et al. (2008); Sobral et al. (2009, 2012a,b), IR: 
Perez-Gonzalez et al. (2005); Wardlow et al. (2011); Mag- 
nelli et al. (2009, 2012); Goto et al. (2010). Such surveys 
have revealed that the redshift range z =1-3 is an essen- 
tial and intriguing era for the study of galaxy formation and 
evolution, since star formation rate density (SFRD) in the 
Universe gradually increases toward z ~ 3 from z >6, has a 
peak at z ~l-2, and decreases sharply from z ~ 1 toward 
2 ~ (e.g., Madau et al. 1996; Hopkins & Beacom 2006; 
Sobral et al. 2012b). The Ha luminosity is widely used to 
derive the star formation activity of galaxies at z < 2, but 
other emission-lines at the bluer end of the galaxy spectral 
energy distribution, such as the [O II] line, are also suffi- 
ciently bright to be widely used and can be employed up to 
z ~ 4, after being calibrated relative to Ha. The Ha emis- 



© 2012 RAS 



2 M. Hayashi et al. 



sion line is a robust star formation indicator which has been 
well-calibrated with data in the local Universe (e.g., Kenni- 
cutt 1998), and is significantly less affected by dust extinc- 
tion than bluer emission lines, such as [O II] . However, the 
Ha line is redshifted into near-infrared (NIR) wavelengths 
for galaxies at z > 0.4, while the [O II] line can be observed 
with an optical instrument until z ~ 1.7, and thus many 
studies/surveys of star-forming galaxies at z > 0.4 rely on 
[O II] luminosity. Unfortunately, while [O II] luminosity is in 
general correlated with the star formation activity, it is also 
dependent on the metal abundance and the ionization state 
of nebular gas. The indirect relation with the star formation 
activity complicates the estimation of SFR from [O II] lumi- 
nosity. Nevertheless, it is empirically calibrated and exten- 
sively utilised as a important SFR indicator for galaxies at 
2 > 1 (e.g., Kennicutt 1998; Moustakas et al. 2006; Gilbank 
et al. 2010). 

In order to better relate the [O II] luminosity of a galaxy 
to its SFR, one can study the ratio of [O II] to H a luminosi- 
ties in local star-forming galaxies. However, it is not obvious 
that these calibrations can be blindly applied for galaxies at 
higher redshifts. It is therefore important to investigate their 
validity at earlier epochs in the Universe by studying higher 
redshift galaxies directly. According to recent studies, the 
[Oil] /Ha ratio of galaxies up to z ~ 1.0 is on average con- 
sistent with that of local galaxies (Moustakas et al. 2006; 
Weiner et al. 2007). Sobral et al. (2012a) investigated the 
[O Il]/H a ratio for star- forming galaxies at z = 1.47 by using 
wide-field, deep imaging with two narrow-band filters which 
can catch Ha and [Oil] emission from z = 1.47 galaxies 
simultaneously. This study finds relatively little evolution 
in the line ratio when compared to lower redshift. However, 
the detection of both emission lines for galaxies at z ~ 1.5 
is limited to relatively luminous galaxies with Ha luminos- 
ity larger than ~ 10 42 erg s _1 (i.e., SFR > 10 M yr _1 ), 
and little is known about the line ratio and nature of the 
more numerous fainter star-forming galaxies. Fortunately, 
the [O II] data obtained with the Subaru telescope can reach 
significantly lower luminosities in [O II] and, using a stacking 
analysis, the [Oil] /Ha line ratio can now be studied down 
to galaxies with significantly lower SFRs. 

The structure of this paper is as follows. The data used 
in this paper are described in § 2. The stacking analysis pro- 
cedure and the results from the stacked images are shown in 
§ 3. In § 4, we show the ratios of H a to [O II] luminosities 
for galaxies at z = 1.47, and compare them with those of lo- 
cal galaxies. We then discuss their implications with respect 
to the evolution in dust extinction and metallicity of [O Ill- 
selected galaxies at z = 1.47. Our conclusions are presented 
in § 5. Throughout this paper, magnitudes are presented in 
the AB system, and we adopt cosmological parameters of 
h = 0.7, Q m = 0.3 and fi A = 0.7. 



2 OBSERVATIONS AND SELECTION 

The bulk of the data used in this paper have been obtained 
as part of the High Redshift Emission Line Survey (HiZELS) 
targeting Ha (Geach et al. 2008; Sobral et al. 2009; Best 
et al. 2010), and as part of a matched [O II] follow-up survey 
(Sobral et al. 2012a, b). The details of the data are described 
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Figure 1. The upper panel shows the response curves of the NBh 
(red solid line) and NB921 (blue dotted line) narrow-band filters 
as a function of redshift which corresponds to the wavelengths 
where H a and [O II] lines are detectable, respectively. The lower 
panel shows the bias of observed line ratio of Ha/[0 II] as a func- 
tion of redshift, which is caused by the slight differences of filter 
profiles. The grey region shows the lc dispersion. Since the pro- 
file of the NB921 filter is slightly broader, the Ha/[0 II] ratio can 
be underestimated for galaxies at the edge of the redshift range 
where the NB filters can catch the emission lines. However, since 
most of the galaxies selected with the NB filters are expected to 
be located at redshifts near the peak of filter profile, it is likely 
that the difference of filter profiles, which we correct for by 3%, 
does not significantly influence the measurement of the line ratio. 

in two papers by Sobral et al. (2012a,b), and therefore we 
only summarise them briefly here. 

2.1 Ha and [O II] emitters at z = 1.47 

We use two narrow-band imaging datasets to obtain our 
samples of galaxies at z = 1.47 in two distinct square de- 
gree areas: the UKIDSS Ultra Deep Survey (UDS: Lawrence 
et al. 2007) and the Cosmological Evolution Survey (COS- 
MOS: Koekemoer et al. 2007; Scoville et al. 2007) fields. 
NBh narrow-band data (A c = 1.617(i.m and AA = 0.021u.m) 
were taken with the wide-field camera (WFCAM; Casali 
et al. 2007) on UKIRT, while NB921 narrow-band data 
(A c = 9196A and AA = 132A) were taken with Subaru 
prime focus camera (Suprime-Cam; Miyazaki et al. 2002) 
on the Subaru Telescope (Figure 1). Combining with broad- 
band imaging data in H and z' which cover the same wave- 
length range as the individual narrow-bands, to estimate the 
continuum level of the spectrum underlying the emission 
line, the narrow-band imaging can measure H a and [O II] 
emission, simultaneously, for galaxies at z = 1.47 (Sobral 
et al. 2012a,b). 

The NBh and H images in the UDS field cover an ef- 
fective area of 0.78 deg 2 where the regions with bad qual- 
ity caused by cross-talk and bright stars are masked, while 
the images in the COSMOS field have an effective area of 
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1.6 deg 2 . The 5cr limiting magnitudes in NBh are ~22.1 
and ~21.9 in the UDS and COSMOS fields, respectively, al- 
though the depth is slightly dependent on the position (see 
Figure 2 and Table 2 in Sobral et al. (2012b) for the details). 

The NB921 data in UDS field are drawn from the 
archive of Suprime-Cam (Ouchi et al. 2010) and reduced 
as described in Sobral et al. (2012a), while public data are 
available for z (Furusawa et al. 2008). In the COSMOS field, 
the NB921 data were taken in service mode with the Sub- 
aru telescope in December 2010 (Sobral et al. 2012b), and 
they cover 69% of the region where the NBh and H data 
are available (i.e., 1.1 deg 2 ). The 5<r limiting magnitudes in 
NB921 are 25.8 and 24.0 in the UDS and COSMOS fields, 
respectively. 
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Catalogues of Ha emitters at z — 1.47 in both the UDS 
and COSMOS fields are presented in Sobral et al. (2012b). 
They include 188 Ha emitters in the UDS down to an Ha 
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Ha luminosity log(L<Ha/erg s 

at z = 1.47) and 325 Ha emitters down to 2.7 x 10 
erg s" 1 cm" 2 (log(L HQ /erg s -1 )=41.56) in the COSMOS 
field, respectively. The number densities of the H a emitters 
are 6.7 and 8.2x 10 -2 arcmin -2 above the flux limits in UDS 
and COSMOS fields, respectively. 
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The catalogues of [On] emitters at z = 1.47 are up- 
dated from those of Sobral et al. (2012a). Emission line 
galaxies are selected as galaxies with colour excess larger 
than 5a (i.e., E > 5) and equivalent width (EW) larger 
than 25A, which corresponds to EWo > 10A in rest frame 
at 2 = 1.47, where photometry is conducted with 2" diam- 
eter aperture. Then, [Oil] emitters at z — 1.47 are identi- 
fied based on the colours, photometric redshift, and spectro- 
scopic redshift (c.f. Sobral et al. 2012b). We note that the 
spectroscopic redshift information indicates that a signifi- 
cant number of galaxies at z ~ 1.3, selected as narrow-band 
emitters due to strong Balmer/4000A break, are incorrectly 
included as [0 11] emitters. However, most of these contami- 
nants are removed from the sample by filtering out galaxies 
with EWo < 20A ni-z > 0.55. Thus, the numbers of [O II] 
emitters selected are 2735 down to 4.2 x 10 -18 erg s -1 cm -2 
(log(L [0 ii]/erg s -1 )=40.76) in the UDS field, and 718 down 
to 1.3 x 10" 17 erg s" 1 cm -2 (log(L [0 n]/erg s -1 )=41.24) in 
the COSMOS field, respectively. The number densities of 
the [Oil] emitters are 0.97 and 0.18 arcmin -2 above the 
flux limits in the UDS and COSMOS fields, respectively. 

By matching the [Oil] emitters with the Ha emitters, 
Figure 2 shows the relation between [O II] and H a luminosi- 
ties for the [0 11] emitters in both the UDS and COSMOS 
fields, and the fraction of [O II] emitters with a significant 
H a emission detection as a function of [O II] luminosity. As 
described in §1, the dual emitters with detections in both 
[O II] and H a are mainly limited to galaxies with [0 11] lu- 
minosities higher than ~ 10 42 erg s -1 . For [Oil] emitters 
with [Oil] luminosities lower than ~ 10 42 erg s -1 , the frac- 
tion is less than ~ 10%, clearly showing that we can only 
investigate the relation between [O II] and H a luminosities 
individually for a small fraction of faint [O n] emitters. 



Figure 2. (Upper panel) The variation of Ha luminosity as a 
function of [Oil] luminosity for the [Oll]+Ha dual emitters at 
z = 1.47, including the upper limits of the Ha luminosity for the 
[O II] emitters without a detection of H a emission. Red symbols 
show the [O II] emitters in the UDS field, and blue ones show 
those in the COSMOS field. The lines show the ratio of H a to 
[Oil] for different levels of dust attenuation (Sobral et al. 2012a). 
(Lower panel) The fraction of [O II] emitters with H a emission de- 
tected as a function of [O II] luminosity. Note that the COSMOS 
fractions are typically higher since the NBh data are Ri 0.2 dex 
deeper in that field. The arrows indicate an upper limit to the 
fraction when there is no dual emitter in the luminosity bin. Er- 
rors in the fraction are estimated based on Poissonian statistics. 
Note that the dual emitters with detections in both [O II] and H a 
are mainly limited to galaxies with [O II] luminosities higher than 
~ 10 42 erg s -1 . 

2.2 Differences between the profile of the two 
narrow-band filters 

Although the NBh and NB921 filters are well matched and 
enable us to detect both H a and [O II] emission lines from 
individual star-forming galaxies at z = 1.47, there is a slight 
difference in the redshift coverage (Figure 1). The redshift 
coverage for [O II] lines detected by NB921 filter is a slightly 
wider than that by NBh filter. It is possible that this differ- 
ence causes us to underestimate the H a luminosity for [O II] 
emitters at redshifts corresponding to the edge of NBh fil- 
ter where the transmission is lower than that of NB921. To 
quantitatively evaluate how accurately the H a/[0 II] line ra- 
tio can be measured by the two narrow-band images and how 
this effect may influence the results of the stacking analysis, 
a simulation was conducted following Sobral et al. (2012a) 
(see §4.5 of their paper for more details). In summary, we 
make a sample of galaxies at z =1.40-1.52 with intrinsic line 
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ratios between and 2.0, where the distribution of [O II] lu- 
minosity that galaxies have is based on the [O II] luminosity 
function, and then investigate the measured Ha/ [Oil] line 
ratios with the two NB filters. Figure 1 also shows the re- 
sult of the simulation: the bias (i.e., the difference between 
a measured Ha/ [Oil] and an input Ha/ [Oil] ratio) is dis- 
tributed around unity in most of the range of the filter pro- 
file, although in the edge of the profile the observed line ratio 
is biased towards lower values. However, we make sure that 
the average value of bias is fairly close to unity. Because we 
focus on the average ratio of H a/[0 II] by a stacking analysis 
in this paper, the slight difference of filter profile has only a 
small effect on the results over the whole investigated range 
of [O II] luminosity. To further quantify this, we investigate 
the bias values for galaxies with a given range of observed 
[O II] luminosity similar to luminosity bins shown in Table 
1; we find that the [O Deselected galaxies have a redshift 
distribution peaking at z ~ 1.46 and with most sources hav- 
ing redshifts where the filter response is high, leading to an 
average value of the bias of ~0.97 in all bins, namely close 
to unity again. Thus, although the Ha/[Oll] ratio can be 
underestimated by 3% (which we correct for), we conclude 
that the difference of the profiles in the two NB filters gives 
no significant influence on the results we find in this paper. 



2.3 Stellar mass for emitters at z = 1.47 

Stellar masses are estimated by an SED-fitting method for 
Ha and [Oil] emitters at z = 1.47 in both fields consis- 
tently (Sobral et al. in preparation). The details of the proce- 
dure in the SED-fitting are described in Sobral et al. (2011). 
The SED templates are created with the stellar population 
synthesis model by Bruzual & Chariot (2003) and Bruzual 
(2007) under the assumption of a Chabrier (2003) initial 
mass function (IMF) and a range of exponentially declin- 
ing star formation histories. In addition, the dust extinction 
law of Calzetti et al. (2000) is used in the SED-fitting. Note 
that we have confirmed that the mass estimates are robust 
against the contribution of emission lines to the broad-band 
photometry, as they are based to fits to ~ 30 photomet- 
ric bands and hence line contamination in one or two bands 
does not have a strong influence. Among [O II] emitters, stel- 
lar masses are derived for 2708 (99.0% of the sample) and 
706 (98.3%) [O II] emitters in the UDS and COSMOS field, 
respectively. The others, i.e., ~ 1% of the samples, do not 
have reliable stellar mass due to large errors in photometry 
and/or non-detection in most bands, suggesting that such 
[O II] emitters are likely to be less massive galaxies. Galax- 
ies for which the SED fitting fails are found to be very faint 
in the rest-frame optical (i.e., AT-band), which supports the 
indication that they are likely to be low-mass galaxies. In 
this paper, we restrict analysis to [Oil] emitters with stel- 
lar masses larger than log(M*/M0)=9.5. This mass cut is 
made in order to maximise the number of [O II] emitters in- 
cluded in the sample, but without picking up many lower 
mass galaxies for which the survey becomes incomplete. In 
addition, the SDSS sample described below, which is used for 
the comparison of galaxies at z = 1.47 with local galaxies, is 
also highly complete above this mass limit. After applying 
the mass cut, the numbers of [O II] emitters are 643 and 212 
in the UDS and COSMOS fields, respectively. 
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Figure 3. The left panel shows stellar mass as a function of 
[Oil] luminosity. Red dots show [Oil] emitters at z = 1.47 in 
the combined sample and blue dots show galaxies in the SDSS 
sample with similar median mass distribution to that of [O II] 
emitters at z = 1.47. Magenta circles and cyan diamonds show 
the median stellar mass of galaxies in each [O II] luminosity bin 
for the samples of [Oil] emitters and SDSS galaxies, respectively. 
The right panel shows mass distribution of each samples: the black 
thin histogram shows the whole sample of SDSS galaxies with 
log(L[ OII ]/erg s~ 1 )>41 and EWq > 10A at z ~ 0.1, the blue 
thick one shows those of local galaxies extracted from the SDSS 
sample, and the red thick histogram shows the mass distribution 
of [Oil] emitters at z = 1.47 in the combined sample including 
the UDS and COSMOS fields. 



2.4 A comparison sample of galaxies at z ~ 0.1 

We also use SDSS DR7 data to compare the results obtained 
for [Oil] emitters at z = 1.47 with local galaxies (Abazajian 
et al. 2009). Galaxies in the redshift range of z =0.07-0.1 are 
extracted from the SDSS spectroscopic catalogue. This red- 
shift range is chosen so that the galaxies have small enough 
apparent angular sizes that most of their light is included 
in the fibre whilst guaranteeing that the sensitivity is still 
very high. An aperture correction on the emission line lu- 
minosities is still required, since the spectroscopic measure- 
ments are done with 3"-diameter fibres. The fractional flux 
loss from the fibre is estimated from the ratio of total mass 
to fibre mass according to the same procedure adopted in 
Sobral et al. (2012a). Moreover, note that the SDSS spec- 
trophotometric calibration takes account of any wavelength 
dependence of the seeing, so the emission line flux ratios are 
also unaffected by this (Adelman- McCarthy et al. 2008). 

By applying our selection criteria, we have selected 8285 
SDSS galaxies with log(L [0 n]/erg s _1 )>41 and EW > lOA. 
Among these, only galaxies with detected H a emission are 
used. However, we note that all but 38 (0.46%) galaxies 
have H a line detected, implying that the removal of galaxies 
without H a detection gives no significant bias in the study 
of Ha/[Oll] line ratios. We are aware that the SDSS spec- 
troscopic sample is magnitude-limited, i.e. roughly mass- 
limited. However, as described in §2.3, the criterion of 
log(M*/M Q )>9.5 is applied to both the SDSS sample and 
the [O II] emitter samples at z — 1.47 so that high complete- 
ness is kept above the mass limit for both samples; this was 
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another factor driving our choice of redshift range for the 
SDSS sample. After applying our mass limit we obtain our 
SDSS sample, which contains 7271 galaxies. 

The SDSS sample is used to compare the Ha/ [Oil] ra- 
tios of [O II] emitters at z = 1.47 with those of local galaxies 
at 2=0.07-0.1 in §§3 and 4. We therefore must assure that 
both samples are fully comparable and that differences be- 
tween them are not arising from a different distribution in 
e.g. stellar mass. Since correlations between SFR and stellar 
mass have been claimed at z ~ 2 (e.g., Daddi et al. 2007), 
it is possible that [Oil] luminosity as well as Ha luminos- 
ity and metallicity are also correlated with stellar mass for 
the z = 1.47 [Oil] emitter samples. We thus use the SDSS 
sample to construct a sample of local galaxies with a similar 
distribution of mass to that of [Oil] emitters at z = 1.47 in 
order to reduce any mass-dependent bias on this study and 
so allow a robust comparison. Galaxies in the SDSS sample 
are selected at random so that the mass distribution is the 
same as that of [On] emitters at z = 1.47 in each bin of 
[O II] luminosity as shown by Figure 3. The SDSS sample 
tends to include galaxies with fainter [O II] luminosity, while 
[Oil] emitters at z = 1.47 have brighter luminosity. Since 
the mass distribution is dependent on the [O II] luminosity, 
the histograms of stellar mass for all galaxies in each sample 
are not in perfect agreement. However, the median masses 
in each luminosity bin are in good agreement (Figure 3). 
In the following sections, we use the matched SDSS sample 
containing 1656 galaxies at 2=0.07-0.1 to compare with the 
results for [Oil] emitters at z = 1.47. 



3 ANALYSIS AND RESULTS 

In order to investigate the relation between Ha and [On] 
luminosities for galaxies at z = 1.47, we stack and average 
the NBh images from which the corresponding H image is 
subtracted (clipping the pixels which deviate more than 3a), 
and then estimated the average H a flux for [O II] emitters at 
z = 1.47 in the UDS and COSMOS fields. Note the results 
do not change if a median stacking is carried out instead. 



3.1 Stacking method 

For all the samples of galaxies to be stacked, a 80" x 80" re- 
gion around individual galaxies is extracted from both NBh 
and H images. Some galaxies are rejected from each sam- 
ple, because they are found in masked regions in the NBh 
image, where some artefacts and regions with low quality 
caused by bright/saturated stars are seen. Then, the zero 
point of magnitude and the point spread function (PSF) 
are matched between the NBh and H images for individ- 
ual galaxies. The original NBh and H images are scaled to 
have same zero point of magnitude. We note that Sobral 
et al. (2012b) show that small colour- correct ions in (H— 
NBh) colour are required to trace the underlying contin- 
uum level and estimate emission line flux correctly, simply 
because the NBh filter is not in the centre of the if-band. 
We thus performed the correction using the (J — H) colour 
according to Sobral et al. (2012b). The PSF in the original 
NB H image has a FWHM of ~ 0.8", while that in the H 
image has FWHM of ~ 1.0". The NB H PSF is therefore 
degraded to match the iif-band seeing by the convolution 
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Figure 4. Comparison between H a luminosities measured on the 
stacked (NBn-ff) image and median values of Ha luminosities 
extracted from the catalogues, for H a emitters which are individ- 
ually detected as galaxies with (//-NBh) colour excess. We plot 
both the individual fields and our combined sample. The median 
and standard deviation of AlogLna, the difference between both 
Ho luminosities, arc 0.018 and 0.045, implying good agreement 
between the two H a luminosities, which confirms the validity of 
our stacking analysis procedure. 



of a Gaussian kernel. After performing these corrections for 
the zero point of magnitude and PSF, (NBH-i?) images are 
created by the subtraction of the H image from the NBh 
image. 

The individual (NBh-H) images are averaged together 
to make the stacked (NBh-^/) image for each sample. The 
photometry is calculated at the centre of the stacked im- 
age with various apertures ranging from 2.0" up to 6.0" 
diameter. We note that the photometry with 4.8" diameter 
aperture is found to recover the total flux of the emission 
line on the stacked image, and then the (NBh-H) colour is 
converted into an emission line flux. Tuned aperture sizes 
could in principle be dependent on the emission line flux, 
but we do not find the significant dependence between the 
size and the luminosity. Thus, we use a common aperture 
with 4.8" diameter to measure the total fluxes in all sam- 
ples. We also note that both Ha and [Nil] emission lines 
can contribute to this emission flux, but no information of 
[Nil] flux for individual emitters is available. Therefore, we 
corrected for the contribution of [N II] flux by assuming that 
[Nll]/Ha ratio of the emitters is 0.22 (Sobral et al. 2012a). 

The [N Il]/H a ratio is known as a indicator of metallic- 
ity, and thus it is probably dependent on the stellar mass. 
Indeed, we find for the SDSS galaxies that the [N Il]/H a line 
ratio correlates with stellar mass more strongly than with 
H a or [O II] luminosities. Our samples include the emitters 
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with a wide range of stellar mass, but the median stellar 
mass in each sample only varies by a factor of four, from 
5 x 1O 9 M to 2 x 1O 1O M . According to the mass-metallicity 
relation that Yabe et al. (2012) have found for star-forming 
galaxies at z ~ 1.4, these median stellar masses correspond 
to the [Nil] /Ha ratios of 0.14-0.25, which are comparable 
to or slightly lower than the value we assumed. Even if the 
lowest ratio of [N Il]/H a=0.14 is used instead, the H a lumi- 
nosities would be increased by only a factor of 1.07. It should 
be noted that the assumption of constant typical value of 
[Nil] /Ha is reasonable, since average Ha luminosities are 
derived in this study (see also Figure 4). 

Contamination in the [O II] emitter samples would lead 
to the Ha fluxes measured on the stacked images being 
underestimated, since the images without Ha flux for the 
contaminants are also stacked. Although the bulk of the 
contaminants are removed using spectroscopic and photo- 
metric redshifts and colour cuts (§2.1), some contamination 
will remain. Using the spectroscopically confirmed NB921 
emitters, and accounting as best as possible for biases in 
the selection of spectroscopic targets, the residual rate of 
the contamination is estimated to be around 15±7%. More 
specifically, it is found that the z = 1.47 [O II] emitter sample 
could include ~ 5% contaminations at z < 1.0 and ~ 10% 
ones at z > 1.0 (most of which are galaxies with a strong 
Balmer/4000A break at z ~ 1.3). We have also found that 
there is no clear trend of the contamination rate with emis- 
sion line flux. Thus, a correction for our best-estimate of the 
contamination, 15%, is applied to the stacked fluxes for the 
[O II] emitter samples in an attempt to account for this. 

AGNs may be included in our emitter samples as con- 
taminants. However, galaxies hosting a strong type-1 AGN, 
with their SED dominated by AGN light, are likely to fail in 
our SED fitting procedure because our code does not include 
an AGN template SED; they will thus be excluded from our 
analysis as only galaxies well-fitted by galaxy SED templates 
are used. Galaxies hosting an obscured type-2 AGN may be 
included as they can be fitted with galaxy SED templates. 
However, such type-2 AGNs usually have a Ha/[Oll] ratio 
above unity (e.g., Cid Fernandes et al. 2010), and so any 
AGN contamination in the z=1.47 [Oil] samples would, if 
anything, bias the stacked Ha flux upwards, making the 
underlying trend found in Section 3.2 even stronger. 

To estimate the error in the flux measured on the 
stacked image, 10,000 identical apertures are distributed at 
random across blank sky region around the galaxy, and then 
the la error is derived by fitting a Gaussian profile to the 
histogram of the sky counts. 



3.2 Stacking analysis for the emitter samples 

Samples of galaxies to be stacked are created from the Ha 
or [O II] emitters in the UDS and COSMOS fields by divid- 
ing them on the basis of their emission line luminosities, so 
that they contain nearly equal numbers of galaxies. How- 
ever, the number of galaxies for the sample in the brightest 
(faintest) luminosity bin is smaller (larger) than the other 
samples. The stacking analysis was conducted firstly for Ha 
emitter samples. The analysis for galaxies with H a emission 
detected individually enables us to check our stacking analy- 
sis by comparing our results with individually measured H a 
luminosities. Figure 4 shows the comparison between median 
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Figure 5. The comparison between Ha and [Oil] luminosities 
for [O II] emitters at z = 1.47. We show the mean H a luminosities 
measured on the stacked (NBjj-ff ) images as a function of median 
[Oil] luminosities in the UDS and COSMOS fields. We also plot 
the combined sample. The gray scale map shows the distribution 
of SDSS galaxies at 2=0.07-0.1, and blue filled circles are median 
H a luminosities for the local galaxies in each [O II] luminosity 
bin. In this plot, the Ho luminosities are corrected for all effects 
discussed in §§2.2 and 3.1. 



Table 1. The average [Oil] and Ha luminosities for [Oil] emit- 
ters in the combined sample in the UDS and COSMOS fields. 
The first and second columns show the median [O II] luminosity 
and the number of [O II] emitters to be stacked in each sample. 
The third and fourth columns show the mean Ha luminosities 
measured on the stacked (NBjj-ff) images and the flux ratios of 
Ha to [Oil]. The contribution of [Nil] to the flux measured on 
the stacked image is removed under the assumption that the typ- 
ical [Nll]/Ha ratio is 0.22 (Sobral et al. 2012a). The flux ratios 
are corrected for any effects discussed in §§2.2 and 3.1. The [Oil] 
luminosities which are originally measured with a 2" aperture are 
corrected by a factor of 1.6 to convert to the total luminosities 
(§3.2). 



[O II] luminosity 


Number H a luminosity 


flux ratio 


(erg s- 1 ) 


(erg s- 1 ) 


(Ha/[Oll]) 



42.22^0 \l 
42 01+ ° :os 
41 90+ 06 

+0.07 
' -0.08 
41 50 +0 19 

4U.DU_ 2g 



61 
127 
125 
126 
124 
246 



42.40 



+0.01 
0.01 
+ 0.01 



42 1 V 
4Z - id -o.oi 

42 02 +a02 
^■ u ^_0 .02 

41 96+ 02 

41 65 +a04 
-8.Q5 
41 53+ - 04 



1.15±0.35 
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0.92±0.13 
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H a luminosity drawn from the catalogues for each H a emit- 
ter sample and that measured on the stacked image. The me- 
dian and standard deviation of Alog(LH Q /erg s _1 ), the dif- 
ference between both Ha luminosities, are 0.018 and 0.045, 
respectively. The good agreement shows that our procedure 
in the stacking analysis recovers well the average of actual 
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individual measurements. We then moved on the stacking 
analysis for samples of [O II] emitters which are classified on 
the basis of their [O II] luminosities. 

Figure 5 shows the average H a luminosities of galax- 
ies at z = 1.47 in the individual fields of UDS, COSMOS 
and the combination of the two, as a function of their [O II] 
luminosity (see also Table 1). The Ha luminosity is esti- 
mated from the stacked NBh - H image with a 4.8"aperture 
(§3.1), and the [Oil] luminosity is a median value of indi- 
vidual [O II] luminosities of emitters in the sample which are 
derived from NB921 imaging. It should be noted that we ap- 
ply an aperture correction to the [O II] luminosities in order 
to ensure a matched comparison with the H a luminosities. 
We obtain this correction by comparing the 2" photometry 
(Sobral et al. 2012a,b) and that obtained with 4.8"apertures 
on the stacked NB921-z' images. We find that a correction 
of 1.6 is required to recover total stacked luminosities. Also, 
the H a luminosities are corrected by 3% for the bias caused 
by the difference of profile of the two narrow-band filters 
(§2.2) and by 15% for the possible contamination of emitters 
at different redshifts that might be included in the sample 
(§3.1). 

There is a clear positive correlation between the mean 
[Oil] and Ha luminosities, with almost constant luminosity 
ratio for [Oil] emitters at z = 1.47. This fact suggests that 
the [O II] luminosity can be easily calibrated as a SFR indi- 
cator for galaxies with SFRs down to 1.4 M Q yr _1 , even at 
z ~ 1.5. The distribution of local galaxies extracted from the 
SDSS sample is also shown in the figures, and the median 
H a luminosities are plotted in each [O II] luminosity bin. For 
z — 1.47 [O II] emitters with a given [O II] luminosity, H a lu- 
minosities are on average lower than those of local galaxies, 
implying that the ratio of H a to [O II] is different between 
local galaxies and [O II] emitters at z = 1.47. To evaluate the 
statistical significance of the difference in H a/[0 II] ratio, we 
compare the Ha/ [Oil] ratios of [Oil] emitters at z = 1.47 
with those of local galaxies, calculate how significantly the 
[O II] emitters at z = 1.47 deviate towards a lower value from 
the local galaxies in each luminosity bin, and then combine 
these in quadrature (using equal weighting) to give a total 
offset significance. We find that there is an overall difference 
in the line ratio which is significant at the 2.8a level. 

The systematic offset is intriguing. It is unlikely that 
the offset is caused by a difference in stellar mass, since the 
mass distributions of the two sample are matched to each 
other (§2.4). Possible systematic effects that might drive the 
offset are discussed in the next sub-section, but are unlikely 
to be able to account for all of the difference, suggesting a 
small but genuine cosmic evolution of the mean H a to [O II] 
ratio. 

3.3 Possibility of underestimation of H a 
luminosity 

One possible cause of bias towards low H a luminosities for 
[Oil] emitters at z = 1.47 with a given [Oil] luminosity is 
contamination in the [O II] emitter samples from emitters 
at different redshifts which are not removed by our spectro- 
scopic redshift, photometric redshift and colour-colour selec- 
tion cuts. As discussed in §3.1, a correction for the contami- 
nation, 15%, has already been applied to the stacked fluxes. 
If the upper limit to the contamination of 22% is the true 



contamination, then the calculated H a luminosities for the 
[Oil] emitters can be also underestimated by ~7%. 

Another possibility is that the contribution of [Nil] 
emission is lower than what we assume in this paper as dis- 
cussed in §3.1. This effect can result in the underestimation 
of Ha luminosity by ~7% at most. Finally, there is a bias 
caused by the difference of profile of the two narrow-band 
filters (§2.2) which can lead to the Ha luminosities being 
underestimated in the wings of the filter profile. Our simu- 
lations have led us to apply a ~3% correction factor for this, 
but these simulations assume a uniform redshift distribution 
for the [O II] emitters across the redshifts probed. If large- 
scale structure causes the [O II] emitters to be clustered at 
the low or high redshift end of the filter coverage, the true 
correction factor could potentially be a few percent higher. 

In order for the Ha luminosities at z = 1.47 to be 
comparable to those of galaxies at z ~ 0.1 for a given [Oil] 
luminosity, all of the effects that might potentially cause an 
underestimation of H a luminosities would have to be at or 
above their maximal combined values of ~ 15%. Thus, it is 
likely that the offset of H a luminosities for galaxies with a 
given [Oil] luminosity is real, and the Ha/ [Oil] ratios of 
[Oil] emitters at z = 1.47 seem to be shifted towards lower 
values than those of the local galaxies to some extent. In 
the next section, we discuss the implication of the different 
Ha/[0 n] ratio in the properties of galaxies. 



4 DISCUSSION 

4.1 Line ratios for the local galaxies 

As discussed by Sobral et al. (2012a), Ha/[Oll] ratios can 
be dependent on both dust extinction and metallicity, be- 
cause the two factors both influence the [O II] luminosity 
more than the H a luminosity. Therefore, before discussing 
the line ratios found in the previous section (Figure 5), the 
SDSS sample with matched mass distribution is used to eval- 
uate how the two factors of dust extinction and metallicity 
relate to the ratio of Ha/ [Oil] and to understand which 
factors are main driver of the difference in the line ratio be- 
tween [Oil] emitters at z = 1.47 and galaxies at z ~ 0.1. 
The Balmer decrement of Ha/H/3 is an indicator of dust 
extinction, while the ratio of [Nll]/Ha, the so-called N2 in- 
dex, offers a rough estimator of metallicity. The Ha/ [Oil] 
ratios for the SDSS sample described in §2.4 are plotted 
as a function of Ha/H/3 and [Nll]/Ha in Figure 6, where 
a polynomial function of the fourth degree is fitted to the 
data. The figure implies that the ratio of H a to [O II] is de- 
pendent on both factors in the sense that it increases with 
higher dust extinction and higher metal abundance. 

It is natural to expect that the H a/[0 II] ratios are sen- 
sitive to dust, since a large amount of dust prevents more 
[O II] emissions at rest-frame 3727A from escaping from the 
star- forming regions in a galaxy than H a emissions at rest- 
frame 6563A. On the other hand, the Ha/[Oll] line ratio 
is also dependent on the metallicity as shown in Figure 
6. According to the metallicities estimated from the mass- 
metallicity relations that recent studies have suggested, the 
figure suggests that the difference in metallicity between 
star-forming galaxies at z ~ 1.4 and z ~ 0.1 could account 
for an offset of A log(H a/[0 ll])~ 0.3dex. 
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Figure 6. (a) The Ha/[Oll] ratios for the SDSS sample as a 
function of the Ha/H /3 ratio. The magenta line shows the poly- 
nomial function of the fourth degree which is best fitted to the 
data, (b) Same as the upper left panel, but as a function of the 
[Nil] /Ha ratio. The red arrow shows the metallicity range in- 
ferred from the mass-metallicity relation of star-forming galaxies 
at z ~ 1.4 (Yabe et al. 2012), while the blue arrow shows that for 
star- forming galaxies at z ~ 0.1 (Tremonti et al. 2004). (c) The 
offset of Ha/ [Oil] ratio from the best fitted polynomial func- 
tion shown in the upper left panel, as a function of [Nll]/Ha, 
i.e, metallicity. The arrows are same ones as shown in the upper 
right panel. The magenta crosses shows the median values in each 
bin of [Nil] /Ha ratio, (d) Same as the lower left panel, but as 
a function of Ha/H/3, i.e, dust extinction, where the best fitted 
function shown in the upper right panel is used. 



In reality, dust attenuation and metallicity are likely 
to be correlated. To investigate the independent effect of 
metallicity on the line ratio, Figure 6 also shows the offsets 
of Ha/ [Oil] ratio from the best fitted polynomial function 
shown in Figure 6, as a function of [Nll]/Ha, i.e, metallic- 
ity. A weak correlation is seen in the sense that the offset 
of A log(H a/[0 II]) gets larger with increasing [Nll]/Ha, 
although the dispersion is large. At fixed dust attenuation, 
the difference in metallicity between star-forming galaxies at 
z ~ 1.4 and z ~ 0.1 suggests an offset of at most ~ O.ldex in 
the mean value of A log(H a/[0 II]). This seems to suggest 
that the metallicity is not the major factor on the offset of 
the Ha/ [On] ratio. The same suggestion is also obtained 
with larger sample of SDSS data by Sobral et al. (2012a). In 
contrast, Figure 6 shows the offsets of Ha/[Oll] ratio from 
the best fitted polynomial function shown in Figure 6, as a 
function of Ha/H/?. This figure suggests that the offset is 
larger with higher dust extinction and the dependence of the 
offset on the dust extinction seems to be stronger than the 
dependence seen in Figure 6 for metallicity, although disper- 
sion is still large. This may suggest that dust extinction is 
the more important factor in changing the Ha/[Oll] ratio. 
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Figure 7. Same as Figure 5, but the ratio of Ha to [Oil] as a 
function of [O II] luminosity. The flux ratios are corrected by 3% 
for the bias caused by the difference of profile of the two narrow- 
band filters (§2.2) and by 15% for the possible contamination to 
be included in the sample (§3.1). The dotted lines show the ratio 
of H a to [O II] in the case of each dust attenuation in H a which is 
derived from the best-fitted polynomial function shown in Figure 
6. 



4.2 Implications for dust extinction and 
metallicity 

The average Ha/ [Oil] ratios are plotted in Figure 7 as a 
function of [O II] luminosity. We find that they shift slightly 
towards a lower value at 2.8a significance than those of 
galaxies at z ~ 0.1. According to the results of Figure 6, the 
ratio of Ha/ [Oil] ~ 1 implies that the typical [Oil] emit- 
ters at z — 1.47 are not dusty, metal-rich galaxies. Moreover, 
it is unlikely that the metallicities of the [O II] emitters at 
z — 1.47 deviate significantly from the mass-metallicity re- 
lation that Yabe et al. (2012) have found for galaxies at 
z ~ 1.4. That is, it seems that metallicity is not a main 
cause of the difference in the Ha/ [Oil] ratio, although we 
cannot completely rule out the possibility that the observed 
line ratios are in part dependent on the metallicity, as sug- 
gested by Figure 6. We thus discuss the line ratios in terms 
of the dust extinction. 

The lines in Figure 7 show the Ha/[0 II] line ratio un- 
der the assumption that the ratios are changed due to the 
amount of dust which is expressed as an attenuation in H a, 
Aua (the best-fitted polynomial function shown in Figure 
6). In this case, all of the Ha/[Oll] ratios for the [Oil] 
emitters at z = 1.47 are distributed around the line cor- 
responding to Ahci~0.35, suggesting that [Oil] emitters at 
z — 1.47 are likely to be less subject to dust extinction. 
This result is not surprising at all, since [O II] emissions at 
rest-frame 3727A should be sensitive to dust extinction. Al- 
though this is within expectations, it is important to directly 
confirm with nebular emissions that typical [O II] emitters 
are on average a less dusty population over the wide range of 
the luminosity down to 10 41 erg s _1 , i.e, SFR=1.4 M yr _1 
(Kennicutt 1998) at high redshift of z = 1.47. 
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Moreover, Figure 7 suggests that local galaxies selected 
based on [O II] luminosity tend to be more subject to dust 
attenuation compared with [Oil] emitters at z = 1.47. 
However, the amounts of dust attenuation are less than 
A(Ha)=l (which is a typical value for Ha emitters) in al- 
most all bins of [O II] luminosity. On the other hand, Sobral 
et al. (2012a) have found that Ha emitters at z = 1.47 
shows the average Ha/ [Oil] ratios consistent with those of 
the local galaxies with A(H«)~1. Thus, [Oil] emitters are 
likely to have smaller amount of dust extinction on average 
than Ha emitters at both redshifts of z ~ and 1.47, and 
more interestingly, typical [Oil] emitters at z = 1.47 seem 
to be less dusty than those in the local Universe. 

Recent [O II] emission surveys at z > 1 have found that 
there are [O II] emitters with red colours comparable to red 
sequence on the colour-magnitude diagram. However, the 
fact that [Oil] emitters at z ~ 1.5 typically show little dust 
attenuation implies that the red [O II] emitters are likely to 
be passive galaxies with AGN activity in the core of galaxy, 
not dusty starburst galaxies. Indeed, we confirm that mean 
H a/[0 II] line ratios for [O II] emitters with (z — K)>2.2 are 
comparable to those for blue [O II] emitters by conducting 
the stacking analysis in the same manner as described above 
for the samples classified by the colours. These results sup- 
port the similar conclusion on the population of red [Oil] 
emitters found in the galaxy cluster z = 1.46 (Hayashi et al. 
2011) as well as those in lower redshifts (Yan et al. 2006; 
Lemaux et al. 2010; Tanaka et al. 2012). Moreover, the re- 
sult we have found in this paper highlights that surveys of 
star formation activity based on only [O II] emissions might 
result in the underestimation to some extent, because [Oil] 
emission surveys tend to be biased toward less dusty galaxies 
and therefore will miss the most dusty starburst galaxies. 



5 CONCLUSIONS 

We investigate the mean relation between H a and [O II] lu- 
minosities for [Oil] emitters at z = 1.47 in the UDS and 
COSMOS fields using a stacking analysis which enables us 
to examine the Ha luminosity of galaxies at z = 1.47 even 
if the individual galaxies are too faint to detect both H a 
and [O II] emission lines simultaneously. The [O II] emitters 
at z = 1.47 are selected with the NB921 narrow-band data 
taken with Suprime-Cam on Subaru Telescope, while the 
H a luminosities are measured on the stacked NBh narrow- 
band data taken with WFCAM on UKIRT. 

We find that on average there is positive correlation be- 
tween H a and [O II] luminosities for not only bright galax- 
ies but also faint ones with [O II] luminosity down to 10 41 
erg s~\ i.e, SFR=1.4 Moyr^ 1 . The trend that galaxies with 
higher [O II] luminosities have higher H a luminosities is con- 
sistent with that of the local galaxies, suggesting that [On] 
luminosities can be used as an indicator of SFR even at the 
high redshift of z = 1.47. 

However, we have to use the [O II] luminosities with 
caution to estimate SFRs at z = 1.47 based on the relation 
calibrated with local galaxies. This is because [O II] emitters 
at z = 1.47 show observed Ha/[0 II] line ratios correspond- 
ing to Aua^0.35 and are less subject to dust attenuation 
than the local galaxies selected based on [O II] luminosity. 
Therefore, [O Il]-selected emitters at z = 1.47 are biased 



toward less dusty populations. The use of dust-correction 
relations derived with H a emitter samples may cause us to 
overestimate the amount of dust extinction and hence the 
dust-corrected SFR. At the same time, surveys of star forma- 
tion activity based on [O II] emissions may miss populations 
of dusty starburst galaxies. 

On the other hand, we note a caveat to our interpreta- 
tion of the results in terms of dust extinction only, because 
the Ha/ [On] line ratio is also dependent on the metallicity. 
Hence the discrepancy of the line ratio between [O II] emit- 
ters at z = 1.47 and local galaxies z ~ 0.1 may be explained 
in terms of metallicity difference. Therefore, the possibility 
that the low H a/[0 II] ratio is not only due to the lower dust 
extinction, but also the lower metallicities of [O II] emitters 
at z = 1.47 than local galaxies, cannot be completely ruled 
out. To distinguish the two factors of dust extinction and 
metallicity completely, we require deep spectroscopy to ob- 
tain the nebular emissions from the individual or stacked 
spectra. 
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